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1. Introduction 

The amplitude and phase of a harmònic oscillator are simple and well understood notions 
in classical physics. The corresponding canonical pair of variables is given by the action 
/ and angle <fi satisfying the classical Poisson bracket {/, 0} = 1. Their quantization 
[H El E] turned out to be a very difficult problem for which a satisfactory solution is 
still lacking. Associating - as suggested by the correspondence principle - the number 
operator N with the action variable, it turns out to be impossible to define a canonically 
conjugate Hermitian phase operator <È> such that the commutation relation [N, Ò]_ — il 
holds. As a consequence, many different and more or less satisfactory approaches to the 
quantum phase problem have been put forward. For a review of the various formalisms 
we refer to the literature 

One of the most 'complete and reasonably consistent' theories is supplied by the 
Susskind-Glogower formalism which has 'become something of a Standard formalism for 
discussing phase problems' (citation from [6 ). The main idea is to consider Hermitian 
cosine and sine operators instead of the qüestionable phase operator itself. The idea 
was suggested by Louisell [TÜ] . implemented independently by Susskind and Glogower 
[TT] and generalized by Lerner [r2| IT3] . 

Here we present a generalization of the Susskind-Glogower formalism from a 
different point of view. The Hermitian cosine C^g and sine Ssg operators introduced by 
Susskind and Glogower in terms of the exponential operators, .Ésg an d E$q, are given 
in the Fock basisj by 

1 1 °° 

C SG = -(E SG + 4 G ) = -^(| n )(n + l| + |n+l}(n|), (1) 

71=0 

1 1 - 

S SG = -(È SG - Èl G ) = - J2(\n)(n + l\-\n + l)(n\). (2) 

n=0 

The corresponding eigenvalue equations Csg|c)sg = c|c)sg an d £sg| s )sg = s|s)sg are 
solved, respectively, by the eigenstates jÜ] 

\ c ) SG = J-^T^J2u n (c)\n), \s) SG = J-^ï^£>[/ n ( S )|n), (3) 

n=0 n=0 

where U n (x) for x = c or s are the Chebyshev polynomials of the second kind. According 
to Q the nth Fock component of the eigenstates, (n|x)sG; is proportional to U n (x), with 
an additional factor i n for the sine states. It is in this sense that we say the Susskind- 
Glogower cosine and sine operators are related with the Chebyshev polynomials of the 
second kind. In this paper we extend this relationship to arbitrary sets of orthogonal 
polynomials on the intervall [-1,1]. The restriction to this intervall follows from the 
envisaged interpretation of x as a cosine c or sine s variable. 

| The Fock basis is formed by the eigenstates |n), n = 0, 1,2, . . ., of the number operator N. The 
states obey the relations of orthonormality (n\m) = S nm and completeness X)^=o l n )( n l = 
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In section |2]we summarize the relevant properties of the orthogonal polynomial sets 
on the intervall [-1,1]. In section|3]we introduce the cosine and sine operators related with 
such polynomial sets and discuss their properties. The corresponding cosine and sine 
probability distributions are considered in section HJ The inverse arccosine and arcsine 
operators and the corresponding cosine-phase and sine-phase probability distributions 
are studied in section EJ Numerical results and graphical representations are given in 
sections 0] to El for particular classical orthogonal polynomial sets and for chosen Fock 
and coherent states. We conclude with a summary in section El 



2. Orthogonal polynomial sets — a reminder 

In this section we summarize the relevant properties of the orthogonal polynomial sets 
[T4*l ITÏÏl ITH] on the intervall [-1,1] and give explicit expressions for the corresponding 
classical orthogonal polynomials used as examples. 

Let w(x) be a real-valued non-negative weight function on the intervall x G [-1,1] 
with the property that all its moments exist and are finite 

r+l 

H n = dx w(x) x n < oo, n = 0, 1, . . . (4) 

Let P n (x) denote a real-valued polynomial of degree n. The set or system of polynomials 
{P n (x)}, n — 0, 1, 2, ... is said to be orthogonal with respect to w(x), if 
-i 

dxw(x) P n (x) P m (x) = d n 5 n , m , (5) 

-i 

where d n > is a finite normalization constant and 5 n ^ m is the Kronecker symbol. 
The weight function w(x) determines the set {P n (x)} up to a constant factor in each 
polynomial P n (x). Different sets of orthogonal polynomials are obtained by choosing 
different weight functions. Each set can be obtained from the fundamental set of non- 
negative powers {l,x,x 2 , . . .} by the Gram-Schmidt orthogonalization method. 
The polynomials of a given set satisfy the completeness relation 

and the recurrence relation connecting three consecutive polynomials [Tlj 

P n +i{x) = (a n + x(3 n ) P n (x) - 7 n P n _i(x), n = 0,1,2,... (7) 
with 7q = by definition. The recurrence coefficients are given by 

a dn+i r> ( b n +i b n \ f3 n d n /a \ 

Pn = , OL n = fj n \ 1 , 7„ = , (8) 

a n \0-n+l a n / Pn-l«rt-l 

where a n and b n are, respectively, the coefficients of x n and x n_1 in the polynomial 
Pn(x) = a n x n + b n x n ~ x + ■ ■ ■. Of particular interest are the polynomial sets generated by 
even weight functions, w(x) = w(—x). In this case all moments of odd order vanish and 
the polynomials of even (odd) degree n contain only even (odd) powers of x, implying 
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Table 1. Weight function w(x), normalization constant d n and recurrence coeficients 
/„ and g n for thc Jacobi polynomials P„ (x). 



P^ v) {x) n>-\,u>-\ 



w(x) = (1-xY {l + xf 

2^+ u+1 T{n + ii + l)T{n + v +1) 
n\(2n + fi + v + l)T(n + fi + v + 1) 

4 j (n + ï)(n + \i + l)(n + v + l)(w + \i + v + 1) 

~ 2n + fi + z/ + 2 y (2n + /i + z/ + l)(2n + /i + z/ + 3) 

" {2n + fi + u){2n + fi + u + 2) V° fi + u + 2/ 



b n = and thus a n = for all n. The polynomials then have a defmite parity according 
to the relation P n (— x) = (—l) n P n (x). 

Related with the set of orthogonal polynomials {P n (x)} we now introduce the set 
of functions {p n (x)} (which are not polynomials, in general) defined by 

Pn(x) = y/w(x)P n (x), P n {x) = P n {x)j\fd n , (9) 

with P n (x) the orthonormal polynomials. The functions p n (x) are orthonormal 
-i 

ÚX Pn{x)p m {x) = S ntm , (10) 



-1 

satisfy the completeness relation 

oo 

^2 Pn(x) Pn(x') = 8(X-X') (11) 
n=0 

and the three-term recurrence formula 

^ Pn+l{x) + Pn-l{x) + 9nPn{x) =Xp n (x). (12) 

Here /_i = by definition and the other recurrence coemcients are given by 



fn a n / d n -\-i b n &n+l 



\; , • 0n = • ( 13 ) 

^ Qn+1 V "n °n 

Note that g n = for the polynomial sets generated by even weight functions. Equation 
(H2J) contains only two independent sets of coemcients, {f n } and {g n }, and is actually 
the well known recurrence relation for the orthonormal polynomials P n (x), since the 
overall factor y/w(x) cancels. For the considerations in the next section it is crucial 
that the coemcients of p n+ \(x) and p n _i(x) turn out to be, respectively, /„ and / n _i 
taken from the same set {f n }- 
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The class of orthogonal polynomial sets on the intervall [-1,1] is very large. For 
illustrative purposes we will restrict ourselves to the remarkable sets of the classical 
orthogonal polynomials. The most general classical polynomials on [-1,1] are the Jacobi 
polynomials Pn l ' u \x) with characteristics given in Table ^ for unequal indices (/i v) 
and in Table El for equal indices (fi = v). Table El contains also the entries for the 
Gegenbauer (A = fi + |), Legendre (À = \) and Chebyshev polynomials of the first 
(À — > 0) and second (À = 1) kind denoted, respectively, by Cn (x), P n {x), T n (x) and 
U n (x). For real-valued polynomials the parameters /i, v and A must be real. The 
restrictions // > — 1, v > — 1 and À > — \ follow from the condition of finite moments 
(0J) . The weight functions and polynomials in Table El have a definite parity implying 
Pn{— x) = (— l) n Pn(x), whereas those in Table[T]lead to pn (— x) = (—l) n pn (x). 

Besides the set of functions {p n {x)} we need two other sets related to it by a 
transformation of the variable x. Since x G [-1,1], we may write x = cos(9 c , with angle 
9 C G [0,7r] and (positive) Jacobian siné'c. The functions defined by 

c n (9 c ) = a/ sin 9 C p n (cos6 c ) = A/sin 9 c w(cos9 c ) P n (cos9 c ) (14) 

then satisfy the following orthonormality and completeness relations 



o 

oo 



d9 c c n (9 c ) c m (9 c ) = 5 n>mi (15) 
J2 cn(9 c ) c n (9' c ) = m-6'c)- (16) 

n=0 

Similarly, the transformation x = sin# s , with angle 9 S G [— 7r/2,7r/2] and (positive) 
Jacobian cos# s , leads for the functions defined by 



s n (9 s ) = a/cos^ p n (sin 9 S ) = A/cos^ü^(sm^) P n (sin6' s ) (17) 
to the following orthonormality and completeness relations 

-7I-/2 

à9 s s n {9 s ) s m (9 s ) = 6 n>m , (18) 

'-■k/2 
oo 

^2 s n (9 s ) s n (9' s ) = m-0's)- (19) 

n=0 

The functions c n (9 c ) and s n (9 s ) are related by an angular shift 

Cn(9 c ) = [sn(9 s )]e s =z-e c , s n (9 s ) = [c n (9 c )] 9c= ^- ds . (20) 

The sets of functions {p n (x)}, {c n (9 c )} and {s n (9 s )} provide an orthonormal basis in 
the Hilbert spaces L 2 (— 1,1), L 2 (0, ir) and L 2 (— n/2, vr/2), respectively. 



3. Generalized and extended cosine and sine operators 

In this section we introduce cosine and sine operators related with arbitrary orthogonal 
polynomial sets on the interval [-1,1], discuss their eigenstates and consider various 
operator relations and expectation vàlues. 
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Table 2. Weight function w(x), normalization constant d n and recurrence coeficient 
/„ for the Jacobi P„ (x), Gegenbauer C„ (x), Legendre P n {x) and Chebyshev 
polynomials of the first T n (x) and second U n (x) kind. 

Polynom. w(x) d n f n 



Pn 



p M {l _ x 2y 2^+ 1 [I> + / , + l)] 2 / ( w + l)(n + 2jz+l 



CÍ A) íl-a: 2 ) A - 



2 



n\(2n + 2/i + l)r(n + 2/x + 1) V (n + fi + \){n + fi + 

7r2 1 - 2A r(n + 2A) / (n + l)(n + 2A) 

n!(n + A)[r(A)] 2 ]j (n + A)(n + A + l) 

2 n + l 
2n + 1 



^/(n + i)(n + |) 



T n {l-x 1 )-^ -r n (r = 2,r„>i=l) 



f/n (l-x 2 )l 



7T 

2 



By analogy with the Susskind-Glogower states (jHJ) we define the states 

oo oo 
\ C ) = Pn(c) \U), \S) = ï n Pn(s) \n) , (21) 

n=0 n=0 

in terms of the functions (jüj). They satisfy the relations of orthogonality 

(c| c ') = ó(c-c'), (s\s') = S(s-s'), (22) 

due to the completeness relation (JTTJ, and resolve the identity 

dc|c)(c| = 1, / ds\s)(s\ = 1, (23) 

'-i J-i 

due to the orthonormality relation (fTüj) . The eigenstates of Susskind and Glogower are 

recovered for the Chebyshev polynomials of the second kind, U n (x), in Table El 

We now determine Hermitian operators C (cosine) and S (sine) by requiring the 

eigenvalue equations C |c) = c\c) and S\s) = s \s) to be satisfied on account of the 

recurrence formula ()12j) . We find 



ò= J2{Y( ln){n + ll + ln + 1){nl ) +9nln){n ^ (24) 

n=0 

oo „ 

S = ^{^(|n)(n + l|-|n + l)(n|)+^|7i)(7i|}. (25) 



71=0 
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The operators consist of an (usual) off-diagonal part of lowering È and raising 
operators, and an (unusual) diagonal part E such that 



C 



S = -{È-Èi) + Èo, 



E 



^2 f n \n){n + l\, È ] = ^ fn \n + l)(n\ , È = ^ g n \n){n\ 



(26) 



(27) 



n=0 



n=0 



n=0 



These expressions are a direct consequence of the recurrence relation in orthonormal 
form (|12|). leading to a raising operator that is the Hermitian conjugate of the lowering 
operator. The diagonal part is absent (Eo = 0) for polynomial sets generated by even 
weight functions (g n = 0). We call the operators C and S without diagonal terms 
Í9n = 0) the generalized, and those with diagonal terms (g n ^ 0) the extended cosine and 
sine operators. Similarly, we distinguish generalized and extended exponential operators 
given by CUS and (C — E ) ±i(S — E ) , respectively. The generalized operators may be 
considered as proper gener alizat ions of the Susskind-Glogower operators (corresponding 
to /„ = 1, g n = 0), and as instàncies of the generalized operators introduced by Lerner 
|12[ Hüj . The extended operators, however, define a more general class due to the 
dependence on two coefficient sequences, {/„} and {g n }. The considerations of this 
paper apply to both types of operators, allowing a unified treatment. 

The cosine and sine operators have a simple matrix representation in the Fock basis 
given by the tridiagonal Jacobi matrices 

\ 



C 



( zgo 

fa 



V ■■ 



fo 
2<7i 

h 







h 

292 
h 






h 

2g 3 



( 



S 



2g 

+i/o 





-i/o 
2<7i 

+i/i 






-i/i 

292 
+i/ 2 






-Í/2 

2g 3 



\ 



\ 



(28) 



The matrix eigenproblems (7\l/ (c) (c) = c\I/ (c) (c) [S^( s) (s) = s\P (s) (s)] lead again to the 
recurrence formula ÍI12I) for the components p n (c) [i n p n (s)} of the eigenvector-columns 
*(C)( C ) [^(S)( s )]. Hence, cosine and sine operators can be related also with tridiagonal 
Jacobi matrices of the form (}28|) having as spectrum the interval [-1,1]. Generalized 
(extended) operators then correspond to matrices with vanishing (nonvanishing) entries 
on the main diagonal. 

We now consider various operator relations and quote them for the extended 
operators. The more compact relations for the generalized operators can be obtained 
by setting g n = for all n. We start with the operators squared 



C 2 



oo 

i E^ 2 + /n-i+4^) \n)(n\ 

n=0 
j oo 

+ o f^Sn + 9n+i) (\n)(n + 1| + \n + l)(n\) 



n=0 
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+ 4 E + 2| + |n + 2)(n|), (29) 

n=0 

oo 

S 2 = T^ + tl+^WW 



4 

n=0 



~ 2 5^^ n ^ n (i n )( n+i i - i n + i >( ri i) 

n=0 
^ oo 

- 4 E + 2 I + ' n + 2 >H)- ( 3 °) 

n=0 

It follows that the classical relation C 2 + S 2 = 1 cannot be satisfied for any choice of 
the parameters /„ and g n . In fact, the only choice [T3| which would do this (f 2n = V%, 
Í2n+i — 0, g n = 0) is incompatible with the recurrence formula (JT2J). 
The commutator and anticommutator are given by 



= \ Y,Un-fLi)\n){n\ 

n=0 
^ oo 

_ 2 ^2fn(9n -9n+i) (\n){n+l\ - \n + l)(n\) 

n=0 

oo 

)(|n)(n+l| + |n + l)(n|), (31) 

n=0 

oo . oo 

[C, S] + = 2Y, 9n \n)(n\ - ~ ^ fj n+1 (\n)(n + 2| - \n + 2) <n|) 

n=0 n=0 
j x 

+ 2 f»(9n + 9n+x) (\n)(n + l\ + \n + l)(n\) 

n=0 

oo 

-\Y. f"(9n + 9n+i) (\n)(n + l\-\n+ l)(n|). (32) 

n=0 

Again, it is not possible to obtain the classical relation [C, S}_ = 0. The commutator is 
not even diagonal unless g n = 0. 

Next we consider some (anti)commutators involving the number operator N: 

[N,Ò\-= -i(S-È ), (33) 
[N,S]_ = +i(C-È ), (34) 

oo „ 

fi* ò h = S [i 2n + !) y (l«><« + l· l + i n + 1><«I) + 2n<? n |n)(n|] , (35) 

n=0 

oo „ 

[^>% = E[(2n + l)^(|n)(n + l|-|n + l)H) + 2^ n |n)H]. (36) 

n=0 

The commutators correspond to the classical relations only if Eq = 0. The mode 
equations following from (jïï3*|) and (jïï4*j) with N as 'Hamiltonian' are given by 

Ò+ C = È Q , È+ S = È q . (37) 
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These are the equations of motion for a harmònic oscillator with equilibrium position 
given by E for g n ^ and by the origin for g n — 0. This gives a nice physical and 
distinctive interpretation of the generalized and extended trigonometric operators as 
normal mode and displaced mode coordinates, respectively. This interpretation is also 
supported by the fact that, due to [N, È ]_ = 0, the operators C — E and S — È obey 
the usual commutation relations with N and hence the normal mode equations. The 
operators C and S and their eigenstates are in quadrature, i.e. related by a rotation 
through an angle ir/2 with generator N: 

C = e- 1 ^ S e +i ^, S = e +i ^ C e~^ È , (38) 

\c) = [e-^ls)}.^, \s) = [e +i ^|c)] c=s . (39) 

This explains nicely the factor i n in the nth Fock component of the sine states. Also, 
being a unitary transformation it relates the sine state with eigenvalue s to the cosine 
state with the same eigenvalue c = s instead of c = y/1 — s 2 . This is one reason why 
the operators fail to reproduce certain trigonometric relations. 

We now come to the expectation vàlues for a system described by the density 
operator p. The expectation value of an operator À is given by the trace (A) p = Sp(pÀ), 
its variance squared by a p AA = (À 2 ) p — (À) 2 p and the correlation of two operators À and È 
by a AB = B] + )p—(A)p{B) p . Using the definitions (jHJ) and and the hermiticity 
of p, the expectation vàlues of the cosine and sine operators are obtained as 

oo oo 

9n Pn (40) 

n=0 n=0 

oo oo 

(S)p = ^2fn^(Pn+l,n) +^2gnPn,n, (41) 
n=0 n=0 

where dt(^s) denotes the real (imaginary) part of a complex number. Compared with 
the expressions for g n = 0, the mean vàlues for g n =/ are shifted by a state-dependent 
amount. The expectation vàlues of the operators squared are given by 



-.00 .00 

( Ò ')p = i + fn-1 + ^l) Pn,n + ^ E f-fn+X ^( P n + 2,n) 

n=0 n=0 

00 

+ E f»(9n + 9n+l) ^(Pn+l,n), (42) 
n=0 

^ OO oo 

^ = I + ^-1 + 4 9n) Pn,n ~ ■= /n/n+1 ^(Pn+2,n) 



n=0 n=0 



+ E f n ( 9n + 9n+1 "> ^(Pn+l,n), (43) 
n=0 

and those of the commutator and anticommutator by 

oo 

([C,S]_) p = ^(fn-fn-^P^n 



2 

n=0 
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oo 

- i ^ 9n ~ 9n+1 ^ [^(Pn+l,n) + 3f(Pn+l,n)], (44) 
n=0 

oo oo 
([C,S] + ) p = 2 J^^p n>n + / n / n+ i 9f(p n+2) „) 

n=0 n=0 

oo 

+ 5^ ^ n ^ n + fi ' n + 1 ) l^iPn+l,n) + 9f(Pn+l,n)]- (45) 
n=0 

We conclude the list of expectation vàlues with those involving the number operator: 

oo 

{[N,Ó\-)p= -Í^/n^(Pn + l,n), (46) 
n=0 

oo 

([N,S]_) p = i^/n»(/Wi lB ), (47) 

n=0 

{[N,Ò\+)p = J2[(2n+l)f n U(p n+1 , n )+2ng n p n , n j, (48) 

n=0 

oo 

([N,S] + ) P = ^[(2n+l)f n $f(p n+1>n )+2ng nPntn \. (49) 

n=0 

The expectation vàlues considered so far involve only the neighbouring density matrix 
elements p n>n , p n +i,n an d p n +2,n- This is a consequence of the tridiagonal structure of C 
and 5 and of having considered at most quadratic expressions in these operators. 
Next we consider some uncertainty relations. The cosine-sine relation is 

°cc v P ss>-^[J2 [(-f" ~ fn-i) Pn,n - 2f n (g n - g n+1 ) [U(p n+hn ) + 9f(p„+i,„)] J • (50) 

\n=0 / 

Contrary to the Susskind-Glogower case, the right-hand-side contains not only the 
vacuum element po,o, but all diagonal elements p n<n and, if g n ^ 0, even the off-diagonal 
elements p n +i, n - The number-cosine and number-sine uncertainty relations are 

<?NN a CC>\ W^fn %(Pn+l,n) ] , °NN °SS > \ [Yl ? n ^«+1,") ) ■ ( 51 ) 

\n=0 / \n=0 / 

Here the right-hand-sides do not depend on g n and on the diagonal elements p n>n . 

We now give examples of expectation vàlues for a Fock state \n). The cosine and 
sine operators have equal expectation vàlues 

(n\C\n) = (n\S\n) = g n , (52) 

which vanish for the generalized operators and yield a state-dependent shift (bias) for the 
extended operators. This is the expected result for a normal and a displaced oscillator, 
respectively. The expectation vàlues of the quadratic operators are given by: 

41 = 41 = \(fn + flil 41 = 0, (53) 

(n\[C,S}Jn) = \ (n\[C,S] + \n) = 2gl (54) 




(n\C 2 + S 2 \n) = ^(/ n 2 + /„ 2 _ 1 + 4^), (55) 

(n\[N,C}_\n) = (n\[N,È\_\n) = 0, <r]& = = 0. (56) 

The cosine and sine operators are uncorrelated and their variances are equal and 
independent of g n . The commutators of C and S with have vanishing expectation 
vàlues and the corresponding correlations vanish as well. These reasonable properties 
hold for both generalized and extended operators. If g n — we have in addition 
(n\C 2 + S 2 \n) = 2(Jqq, i.e. the trigonometric relation (n\C 2 + S 2 \n) = 1 is violated 
to the same extent to which the variances squared deviate from their classical value 
\ (see (J75]) below). In the case of the classical polynomials and for large vàlues of n 
(f n — > 1, g n — > 0), the shifts of the mean vàlues approach zero, the variances squared 
and {n\C 2 + S 2 \n) their classical vàlues, and S]_|n) vanishes. This is the expected 

behaviour on the basis of the correspondence principle. 

We show in Figure [U the variances squared tjçç = cr^ for the Fock states |0), 
|1), |2) and |3) as functions of the Gegenbauer parameter A. For À = 1 (Chebyshev 
polynomials of the second kind or Susskind-Glogower case) the variances squared are 
equal to the classical value \ for all Fock states n^O and to j for n = 0. Similarly, for 
À — ► (Chebyshev polynomials of the first kind) they are equal to | for all n ^ 1 and to 
| for n = 1. For À > 1 the variances squared stay below | for all n, decrease to zero for 
very large À and increase with n for fixed À. For n > 2 and À < (0 < À < 1) they stay 
below (above) |, with a maximum value at some point A max G (0, 1) which approaches 
| for large n. Note that vàlues of the variances squared greater (less) | correspond to 
vàlues of (n\C 2 + S 2 \n) greater (less) 1. 

As a next example we consider the coherent states 



a) = e - ^ J2 -7= l n >> (» = l a l e ^)- ( 57 ) 



n=0 
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The expectation vàlues of the cosine and sine operators are now given by 

(a\C\a) = F 1 costp + Gi, (a\S\a) = F 1 eàmp + Gi, (58) 
and their variances squared and correlation by 

4t = (F 2 - Fl) cos 2 tp + ^(F + - F 2 ) + (G+ - 2F 1 G 1 ) cos <p + G 2 - G 2 , (59) 

41 = (F 2 - Fl) sin 2 if + ^(F + - F 2 ) + (G + - 2F 1 G 1 ) siny? + G 2 - G\ , (60) 

a^ s = (F 2 - F 2 ) cosy> siny? + ^(G + - 2F 1 G 1 )(cosyi + siny?) + G 2 - G 2 . (61) 
Here we have introduced the \ct (-dependent functions 

{F i; G ± } = e-H ^{/ n ; ± <7 n+1 )} ' , (62) 

F2 = e " H S /«/»+! , 7= ' =^ > ( 63 ) 

^ n!V(n + l)(n + 2) 

00 i I |2n 

{F ± ; G i; G 2 } = e-l«l 2 £{-(/ 2 ± fl,); g n ; g 2 n } . (64) 

n=0 

Some other expectation vàlues of interest are: 

(a\[C,S]_\a) = i[F_ -G_(cosy? + siny?)], (65) 
(a\C 2 + S 2 \a) = F + + 2G 2 + G + (cosy? + siny?), (66) 
(a\[N,C]Ja) = -iFiSin^, (a|[JV,5]_|a) = iíicosp. (67) 

For the generalized operators all G functions vanish and the relations above take on 
familiar forms. These forms get usually modified for the extended operators by state- 
dependent terms. Thus, the expectation vàlues of C and S get shifted by the same 
amount Gi(|a|). For the classical polynomials and for large vàlues of \a\ we recover the 
classical results (a|G|a) — > cosy? and (a^la) — > siny?, since F\ — > 1 and G\ — > in 
this limit. Figure El displays Fx(|a|) as a function of \a\ for À = — h, ~, 1 and 5. The 
functions increase monotonically with \a\ for A > and develop a bump followed by a 
valley for A < 0. All functions start at zero for |a| = (vacuum state) and approach 1 
(classical limit) for \a\ — > oo; the approach is slower for higher vàlues of A. 



4. Cosine and sine representations and distributions 

In this section we introduce the cosine and sine representations of an arbitrary state and 
discuss the corresponding probability distributions. The cosine and sine distributions are 
considered on their own right as functions of the eigenvalues c and s, without expressing 
the eigenvalues in terms of phases.f The importance of this kind of distributions for the 
quantum measurement of trigo nometric variables has been stressed in |17j . 

f In the dcfinitions (|24|l (J2SJ) of the operators C (S) there is no reference to a phase and labeling the 
eigenvalues c (s) by cos 9 (sin 9) should be considered merely as a parameterization at this stage. 
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Figure 2. i*i(|a|) for À = —4, i, 1 and 5 as function of \a 



The cosine and sine states (J2TJ) can be used to define cosine and sine representations 
of an arbitrary normalized state with Fock components ip n = (n\if)) by 

oo oo 
^ C \c) = (c\ll>) = 5>n(c)Vn, ^ {S) (S) = (é#> = ^(-i)>„( S ) fa. (68) 

n=0 n=0 

The cosine and sine distributions are then given by \tp^ c \c)\ 2 and \ip^\s)\ 2 or, more 
generally, for a mixed state with density matrix elements p m>n by 

oo oo 

V f\ C ) = E PmnPm(c)Pn(c), V® (fi) = E (-i)" 1 "" p m „ p m (s)p„(s) ■ (69) 

m,n=0 m,n=Q 

Since p m (^)Pn(^) is symmetric in m and n, only the symmetric part of the remaining 
factors contribute to the double sums. Taking the hermiticity of the density matrix into 
account, this symmetric part is 3?[p mj „] for x = c and §ft[(— i) m_n p m ,n] for x = s. The 
cosine and sine probability distributions are then given by 

oo oo 
V l C) ( C ) = 5ZPn,nP'(c) + 2 ^lPm,n}Pm{c)Pn{c), 



(70) 



n=0 



m>n=0 
oo 



Pp W (*) = E^P£( a ) + 2 E «lAn.ne^^lMaKW· (71) 

n=0 m>n=0 

The sine distributions follow formally from the cosine distributions by replacing the real 
part of the density matrix elements p mi „ by the real part of the phase-shifted matrix 
elements p m>n exp[— i~(m — n)]. For a large class of states, including the coherent and 
squeezed states, the density matrix elements p m>n have a phase dependence of the form 
exp[i(m — n)(p]. In this case the sine distributions follow from the cosine distributions 
by the substitutions (c, ip) — > (s,<p— |). In other words, the phases 'measured' by the 
cosine and sine distributions are in quadrature as expected. 
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Figure 3. Cosine (and equal sine) distributions VL^ (c) for the Fock states |0), |1) 
and |5), for A = — j and 1, compared with the classical distribution V cX (c). 



The cosine and sine distributions can be used to calculate the expectation vàlues 
of functions F(C) and F(S) according to 

(F(C)) P = J\cF(c)VP(c), (F(S)) P = £\sF(s)Vf\s). (72) 

We now apply the above considerations to the Fock and coherent states. Using 
í)21jl . the cosine and sine representations of a Fock state |n) are given by 

(c\n)= Pn (c), (s\n) = (-i)> n (s), (73) 

i.e. by the same function p n (x) evaluated for x = c and x = s, respectively, with an 
additional factor (—i)™ for (s|n}.| The corresponding cosine and sine distributions are 
equal and given by 

V$(x) = p 2 n (x) = w(x) P 2 n (x)/d n , (x = c, s). (74) 

This is precisely the expression under the normalization integral (JlUj) for n = m, now 
interpreted as the probability distribution to find the value x in the Fock state |n). The 
distribution for the vacuum state is given by the normalized weight function w(x)/fiQ, 
where /i is the zeroth moment (jU); if Pq{x) = 1 then fi = d . The distribution (|74^) 
vanishes at the n zeros of the polynomial P n (x), n > 1, all of which are known fH] to be 

f This is in perfect analogy with the usual coordinate ((q\n) = h n (q)) and momentum ((p\n) = 
(—i) n h n (p)) representations of the nth energy eigenstate of the harmònic oscillator. 
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Figure 4. Cosine (and equal sine) distributions VL· (c) for the vacuum state |0), for 
A=— |, — > 0, \, |, 1 and 2, compared with the classical distribution (c). 

real, distinct and located in the interior of the interval [-1,1]. The distribution therefore 
shows an oscillatory behaviour with n + 1 peaks. It is symmetric for the generalized 
and asymmetric for the extended operators. The behaviour near the end points x = ±1 
depends solely on the value of w(±l), since P n (±l) is finite. In the particular case of 
the Jacobi polynomials the distributions are asymmetric (symmetric) for /i 7^ v (fi = v) . 
The behaviour near the end points is given by (1 — x) M for x — > 1 and by (1 + x) u for 
x — > — 1. Due to the property Pn\—x) = (—l) n pn'^ix) , the distributions with Jacobi 
indices (/z, v) and (z/, /i) are related by reflection. The vacuum state distributions are 
given by w(x)/do, with d from the Tables^and El They are uniform for the Legendre 
polynomial Pq(x) and equal to the classical distribution 

V f\x) = —L·= (x = c,s) (75) 

7TV 1 — X 

for the Chebyshev polynomial Tq(x). The distribution ()75j) corresponds to a uniform 
arccosine or arcsine distribution (see sectionEJ), yields vanishing expectation vàlues and 
variances squared equal to |. 

The general behaviour outlined above can be clearly observed in the figures El to 
where the classical distribution is displayed for comparison. Figure El shows the cosine 
(and equal sine) distributions for the Fock states |0), |1) and |5) and the Gegenbauer 
parameters À = — |, 1. The distributions for À = —\ (À = 1) diverge (vanish) at the end 
points and have n (n + 2) zeros and n + 1 peaks for n > 1 . The quantum distributions 
oscillate around the classical distribution, approaching it for large vàlues of n.\ 

Figure |H shows the cosine (and equal sine) distributions for the vacuum state |0) 
and for different vàlues of À: — |, 0, |, 5, 1 and 2. For A > \ (A < |) the distributions 

| This is reminiscent of the coordinate and momentum distributions for the nth eigenstate of the 
harmònic oscillator, in which case the zeros of the distributions are the roots of the Hermite polynomials. 
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Figure 5. Cosine (and equal sine) distributions V^'(c) for the vacuum state |0) and 
for the Jacobi parameter vàlues (fj,, v) = (— |, — |), (— |, |), (~, — |), (i, |). 



vanish (diverge) at the end points c = ±1 and have a maximum (minimum) increasing 
with À at c = 0. 

Figure shows the cosine (and equal sine) distributions for the vacuum state |0) 
and various Jacobi parameter sets (/x, v). The distributions are asymmetric and show 
the behaviour near the end points as discussed above. The asymmetry of the vacuum 
state distributions is at first glance an unexpected result, but it may be natural for a 
'displaced' oscillator (/i ^ v). By contrast, the distributions in Figure |U corresponding 
to a 'normal' oscillator (fi = v) show the expected symmetry. 

We now come to the coherent states (|57|). Their cosine and sine representat ions are 

2 00 \ n \n .2 00 \ n \n 

(c\a) = y;^=e in > n (c), (s\a) = e"^ V^e^-f'p^s). (76) 

n — 11 v n — II v 



n=0 v ra=0 

The cosine probability distribution is then given by 

°° r lirl 2n °° lirl m+n -i 

W = eH ° |2 E {\pK c ) + 2 E cos [( m " n)íP]Pn,(cK(c)}, (77) 

n=o 1 n! m (>„) Vmlnl J 

and the sine distribution 7^ ,(s) follows therefrom by replacing (c,<p) — > (s, y — f) on 
the right-hand-side. The dependence on the coherent phase is through cos[(m — n)(p] 
(cos[(ra — n)(<p — |)]) for the cosine (sine) distribution. For a fixed phase (p , the states 
with <p — <p and <p = —(p (ep = ip and <p = tt — (p ) then have the same cosine (sine) 
distribution. These are precisely the phases which yield the same value of the cosine 
(c = cos<y9 ) or sine (s = sin<^ ). States having identical cosine (sine) distributions are 
distinguished by their sine (cosine) distributions, so that both distributions are needed. 
Note also that the sine distributions for the states with (p — ^ ± (fo are identical with 
the cosine distribution for the state with ip = ipQ. 

Figure |ü] shows the cosine distributions for À = — |, 1 of the coherent states with 
a = 1, i, 1 + i, (vacuum state). The curves referring to a = i, 1 + i give also the cosine 




c 



Figure 6. Cosine distributions Pf,' (c) for the coherent states with a = 1, i, 1 + i, 0, 
and for A = — 4, 1. Identical distributions for a = 1, — i, 1 — i, 0, respectively. 

distributions for the complex conjugate vàlues a = —i, 1 — i. The sine distributions for 
the same input data are shown in Figure Here the curves for a — 1, 1 +i refer also 
to the sine distributions for a = —1, —1 + i, respectively. The cosine distributions for 
a — 1, i, 1 + i, and the sine distributions for a = i, ±1, ±1 + i, are identical in that 
order. The distributions for A > peak, as expected, around the cosine and sine vàlues 
corresponding to the coherent phase, <p = (c = 1, s = 0), ip = | (c = 0, s = 1) and 
(p = j (c = s ~ 0.71), but develop unnatural dips for A < 0. The latter distributions 
should be considered therefore with care. The peaks of the distributions get sharper 
with increasing vàlues of \a\. 

5. Arccosine and arcsine operators, representations and distributions 

In this section we consider the inverse trigonometric operators and use their 
eigenstates to define arccosine and arcsine representations and corresponding probability 
distributions. The direct sampling of such distributions in the particular Susskind- 
Glogower case has been considered in [T8] . 

The arccosine and arcsine operators are defined by the power-series expansions 

00 / i \fc 1 

^E^t(7)^ +1 - < 79 > 

fc=0 

They are of the generalized or extended type depending on the type of the corresponding 
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Figure 7. Sine distributions Pfns) for thc coherent states with a = 1, i, 1 + i, 0, 
and for A = — 4, 1. Identical distributions for a = — 1, i, — 1 + 0, respectively. 

trigonometric operator. The series are well defined even for the extended operators, since 
the spectrum of the latter is also confined to the convergence intervall [-1,1]. Combining 
(fT8j) and (f?9|) we obtain arccosC* = |1 — arcsinC and arccosS* = |i — arcsinS 1 , which 
are the quantum versions of the classical relation arccosx = | — arcsinx. 

The operators C and (9 S do not commute, but C commutes with C and S with 
S. It then follows that the states |c) and \s) are also eigenstates of C and 6> s 

6> c |c) = arccos c |c), 6> s |s) = arcsin s \s). (80) 

Writing c = cos# c and s = sin# s , and choosing the principal value branches < 9 C < n 
and — | < 9 S < +| , we define the arccosine and arcsine states by 

oo 

\9 C ) = \/sm8 c [\c)} c=cos e c = ^c n {d c )\n), (81) 

n=0 

oo 

\9 S ) = Vcos9 s [\s)} s=sin e s = ^i n s n (éy |n>, (82) 

n=0 

with c n {9 c ) and s n (6' s ) given in (fl·lj) and (|T7j). The states \9 C ) i\9 s )) are eigenstates of 
C (0 S ) as well as of C (S) with appropriate eigenvalues: 

0c\9 c ) = 9 C \9 C ), S \9 S ) = 9 S \9 S ), (83) 

C\9 C ) = cos9 c \9 c ), S\9 S ) = Ún9 s \9 s ). (84) 

They are orthogonal and yield a resolution of unity according to 

(OM = - 0' c ), (OM = 5(9 S - 8' s ), (85) 
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p+tt/2 

d9 c \9 c )(9 c \ = í, / de s \6 s )(6 s \ = í. (86) 

J-ir/2 



The inverse trigonometric operators and their eigenstates are related by 



6c = e" 1 ^ 



-l-6s 

2 



2 



e~^ N } (87) 



\0c) = [e"^ v \e s ))e s=í -o c , \9 S ) = [e+^ \e c )] dc= ^o s . (88) 

Hence, S (0 C ) is unitarily related to [| 1 — C ] ([| 1 — S ]) and not to C (0 S ) itself. 

We have called the states (J8"ïj) and (|8*2^) the arccosine and arcsine states, respectively. 
The terms cosine-phase and sine-phase states would also be appropriate and are actually 
used for the well-known Susskind-Glogower states 



\9 c )sg = J-£sin[(n + 1)0 c ] \n), \6 S ) SG = J- ^i n sin[(n + 1)(- - 9 S )} \n), (89) 

n=0 n=0 

to which our states reduce in the particular case of the Chebyshev polynomials of 
the second kind by using the representation U n (cos6) = sin[(n + 1) 9}/ sin9. The 
states (|5üj) are usually denoted by Icosó'c) (|sin^ s )) to emphasize the eigenvalues of 
C (S) rather than of C (0 S ). Nevertheless, they are normalized to 5(6 — 0') and the 
measure in the resolution of unity is d9, as in our case. From this point of view our 
notation is more consistent with Dirac's convention. Also, in the present context the 
notations | cos 6^) ( | siné^)) for the states (|HT|) ( (jS2j) ) could give rise to confusion with 
the only reparameterized states [|c)] c=COS c ([|s)] s=sin ej not including the square root of 
the Jacobian. Note that 9 C (9 S ) are not true phases but, strictly speaking, the arccosine 
(arcsine) of c (s) and vary accordingly over a 7r-range instead of a 27r-range characteristic 
for a genuine phase. 

The states \9 C ) (\9 S )) in lj%I|l ( (|B5Jl ) have the same structure as the states \c) (\s)) 
in (f2*T|). with c n {9 c ) (s n (9 s )) replacing p n (c) (p n (s)). These replacements act as a rule 
when switching between appropriate trigonometric and inverse trigonometric quantities. 
Thus, the arccosine and arcsine (or: cosine-phase and sine-phase) representat ions of a 
state \ip) are given by (compare (jüSJl ) 

oo oo 

^\6 C ) = (9 c \íj) = Y,cn(9 c ) fa, ÍJ ie °\6 s ) = (9 S \4>) = Y^H) n s n (9 s )^j n . (90) 

71=0 71=0 

The corresponding arccosine and arcsine (or: cosine-phase and sine-phase) probability 
distributions follow from (J70*|) and (fTTj) by the same replacements 

V^\9 C ) = [PP(c)] Pnic) ^ Cn{6ch Vf-\9.) = [Vf\s)] Pn(s) ^ n{ea) . (91) 

For the expectation vàlues of functions F(0 C ) and F(0 S ) we then have 

»+i 



(F(0 C )) P = Ç à9 c F{9 c )Vf°\9 c ) = j dcF(arccosc)Pf)(c), (92) 
(F(0 s )) p = [ +W/2 d9 s F(9 s )Vf°\9 s )= / +1 d S F(arcsin S )Pf (s). (93) 

J-w/2 J-l 
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Figure 8. Arccosine and arcsine variances squared <Jq q = ffg e for the Fock states 
|0), |1), |2) and |3) as functions of A, compared with the classical value n 2 /12 w 0.82. 

These relations reduce to lf73jl if F depends on C (0 S ) via cos C = C (sin S = S) 
only. The classical uniform distributions correspond to 

Vf c \0c) = ~ on[0,7r], V^\6 S ) = - on (94) 

They yield the expectation vàlues (9 c )d = f and (9 s ) c i = 0, equal variances squared 
a e\e c = a te = 71-2 ! /12 and the classical distribution (fTHJ) for x = cosé'c and x = sin# s . 
The arccosine and arcsine representat ions of a Fock state \n) are given by 

(9 c \n) = c n (6 c ), (9 s \n) = (-i) n s n (6 s ), (95) 

and yield the probability distributions (related by 9 C = |- — 9 S ) 

V\%\9 C ) = c 2 n (9 c ), V§'\0.) = b'M. (96) 

Using (J87j) and the binomial expansion formula, we obtain the following relation between 
the arccosine and arcsine moments 

(-1 O k c \n) = (£)* + £(-!)' ( k ) (\) k ~ l (n\ è\\n). (97) 
i=i 

In particular, for k = 1 and k = 2 we find 

(n|0 c |n) = |- (n| 6> c » = - vr (n| 6>» + (n\ 6>», (98) 

implying equal variances squared, cr^ e = Vqq- For polynomials of defmite parity 
the expectation vàlues (n\0 l s \n) vanish for odd powers /. In this case (n\0 c \n) = ~ 
and (n\0 s \n) = 0. The variances squared are shown in Figure |H1 as functions of the 
Gegenbauer parameter À for the Fock states |0), |1), |2), |3). For fixed À and increasing 
n they approach either from above or from below the classical value 7r 2 /12 = 0.8225. 
The behaviour is similar to the one observed in Figure ^ for Gç* c = a^. 




Figure 9. Arccosine distributions T'L· (0 C ) for thc vacuum state |0), for A = — j, 
— > 0, \, h, 1 and 2, compared with the classical value 1/tt « 0.32. 



The arccosine [arcsine] distributions (|96|) for a Fock state |n) vanish at the n roots of 
the polynomial P n (cos9 c ) [P íl (sin^ s )] and must show therefore an oscillatory behaviour 
with n + 1 peaks. At the end points the distributions vanish, stay finite or tend to 
infinity depending on the behaviour of sm.9 c w{cos9 c ) [cos 6^ u>(sin 6^)] for 9 C — > 0,7r 
[9 S — > 7r/2, — 7r/2]. In the case of the Jacobi polynomials this behaviour is given by 
(9 C ) W , (vr - 9 c ) 2u+1 [(§ - 9 s ) 2ll+ \ (f + e s ) 2v+1 ]. For the Gegenbauer polynomials 
(2/1 + 1 = 2u + 1 = 2A) the distributions vanish (diverge) at both endpoints for À > 1 
(À < 1) and stay finite for À — > (Chebyshev polynomials of the first kind). These 
features can be observed in Figure EU showing the arccosine distributions of the vacuum 
state |0) for A = — |, — > 0, |, |, 1, 2. Here the distribution for À — > (corresponding to 
the polynomial T (x)) coincides with the classical uniform distribution ([94)1 . 

The arccosine (^ c |a) and arcsine (8 s \a) representat ions of the coherent state (|57j) 
can be obtained from (|7üj) by replacing p„(c) and p n (s) by c n {6 c ) and s n (9 s ), respectively. 
The corresponding distributions, V^ C \9 C ) and T'S (^«), are then given by the right- 
hand-side of (J77j) with p n (c) replaced by c n (8 c ) and (p n (c),<^) by (s n (9 s ),ip — |), 
respectively. As a result, the comments made in the paragraph following (|77|) apply 
appropriately. Thus, the coherent states with +(p and — cp (cp and 7r — ip ) have 
the same arccosine (arcsine) distributions. Also, the arcsine distributions for | ± cpo 
and the arccosine distribution for ipo agree if the translation 9 C = | — ^ s is taken into 
account. Neither the trigonometric nor the completely equivalent inverse trigonometric 
distributions can resolve this phase ambiguity. The c (s) distributions 'feel' only the 
vàlues of cos<£> (sinyj), and the 9 C (9 S ) distributions, strictly speaking, only those 
of arccos(cos^) (arcsin(sin ep)). None of the distributions is a true phase-sensitive 
distribution. Actually, they cannot be one by construction. 

We show in Figure ITÜ1 the arccosine distributions for the coherent states |a) with 
a = 1, i, 1 + i, (vacuum state) and for À = — |, 1. The corresponding arcsine 
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distributions are shown in Figure ^2 The distributions should be compared with the 
cosine and sine distributions for the same input data in Figures |ü] and Here the 
distributions develop peaks (dips) for A > (A < 0) too. Their position is close to 
9 C = ip in Figure [Tü] and to 9 S = p in Figure HU where (p — 0, | and | is the phase of 
a. If we expect a distribution to peak around the value corresponding to the coherent 
phase, then the distributions with A < should be disregarded. 

We conclude this section by making two remarks. First, we note the identities 

cos 2 C + sin 2 9 C = 1, cos 2 S + sin 2 S = i, (99) 

showing that in the sense of this trigonometric relation the 'right partner' of C — cos C 
(S = sin & s ) is sin(arccos C) (cos(arcsin S)) and not 5 (C). 

Second, the Hermitian operators C and © s can be used to define two distinct 
noncommuting unitary operators 0] by 

Ü c = exp(i6> c ), Ü s = exp(i6> s ). (100) 

Writing exp(ié>) = cos + i sin and using püjl. we obtain 

C + iS = É +(l + i)È = ^[Ü c + Ü s + ÜÍ-Üt], 

C-iS = & + (l-i)Èo = l[Ü c -U a + Ül + Ul]. (101) 

The combinations C ± iS yield pure shift operators, as in the Susskind-Glogower and 
Lerner cases, only for the generalized operators. In the case of the extended operators 
the displaced combinations (C — È ) ± i (S — È ) should be considered instead. 
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Figure 11. Arcsine distributions (0 S ) for the coherent states with a = 1, i, 1 + i, 
0, and for A = — j, 1. Same distributions for a = — 1, i, — 1 + i, 0, respectively. 

6. Summary 

In this paper we have introduced cosine and sine operators which generalize (extend) in a 
specific way the cosine and sine operators of Susskind and Glogower. Our starting point 
was the observation that the eigenstates of the Susskind-Glogower operators in the Fock 
basis are given by the Chebyshev polynomials of the second kind. We have extended 
this relationship to arbitrary polynomials which are orthogonal on the intervall [-1,1] 
with respect to a weight function. Related with each polynomial set there is a pair of 
cosine and sine operators determined by the recurrence coefficients for the orthonormal 
polynomials and with eigenstates given in terms of the polynomials themselves. The 
relationship implies: i) the eigenvalue equations for the cosine and sine operators are 
satisfied on the basis of the three-term recurrence formula for the polynomials, ii) 
the orthogonality and iii) the completeness relations for the eigenstates follow directly 
from the completeness and orthogonality relations, respectively, for the polynomials. 
Two types of operators emerged naturally, termed generalized or extended operators, 
depending on wether the weight function is symmetric or not. The distinction can be 
interpreted physically in terms of normal mode or displaced mode variables. In the case 
of the classical orthogonal polynomials the Jacobi polynomials with equal (unequal) 
indices give rise to generalized (extended) operators. The generalized operators belong 
to the class of operators considered by Lerner and are proper generalizations of the 
Susskind-Glogower operators. The extended operators define a new class of cosine and 
sine operators which, to our knowledge, has not been considered so far in this context. 
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For both types of cosine and sine operators we have introduced lowering and 
raising exponential operators as well as inverse arccosine and arcsine operators. 
The eigenstates of the trigonometric and inverse trigonometric operators have the 
orthonormal polynomials as Fock components, are orthogonal and yield a resolution 
of unity. This allows them to be used in defming cosine, sine, arccosine and arcsine 
representations and corresponding probability distributions for arbitrary normalized 
states. The phases 'measured' by the cosine and arccosine distributions on the one 
side and the sine and arcsine distributions on the other side are in quadrature, as 
expected. None of the distributions, however, is a true phase distribution. The still 
missing eigenstates of the lowering exponential operators of both types are of a different 
nature and will be considered in a forthcoming paper. 
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